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Abstract: In this paper we generalize the construction of locally boosted black brane space 
time to higher derivative gravities. We consider the Gauss-Bonnet term (with coefficient a') 
as a toy example. We find the solution to the a' corrected Einstein equations to first order in 
the boundary derivative expansion. This allows us to find the a' corrections to the boundary 
stress tensor in the presence of the Gauss-Bonnet term in the bulk action. We therefore obtain 
the ratio of shear viscosity to entropy which agrees with other methods of computation in the 
literature. 
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1. Introduction 

The conjectured duality between thermal gauge theory and gravity in one higher dimensional 
AdS spacetime is a useful tool to extract information about strongly coupled large N gauge 
theory. For example different thermodynamic variables like free energy, entropy etc of large 
N SU(N) gauge theory in (3 + 1) dimensions in large 't Hooft coupling regime can be 
conjectured by calculating the thermodynamic variables of black hole in five dimensional 
AdS spacetime. 

If we consider a black object with translation invariant horizon, for example black D3 
brane geometry, one can also discuss hydrodynamics - long wave length deviation (low fre- 
quency fluctuation) from thermal equilibrium. In addition to the thermodynamic quantities 
the black brane is also characterised by the hydrodynamic parameters like viscosity, diffusion 
constant etc. The black D3 brane geometry with low energy fluctuations (i.e. with hydrody- 
namic behaviour) is dual to some finite temperature gauge theory plasma living on boundary 
with hydrodynamic fluctuations. Therefore studying the hydrodynamic properties of strongly 
coupled gauge theory plasma using the AdS /CFT duality is an interesting subject of current 
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research. The energy momentum tensor of a relativistic viscous conformal fluid is given by 
(up to first order in derivative expansion), 



T^ u = (e +p)u fl u u +pr] txv - (f.f) 

where is fluid 4- velocity with u^u^ = — 1, e is energy, p is pressure and r] is shear viscosity 
coefficient. a^ u is defined in 2.13. Conformal invariance implies that e = 3p. The first attempt 
to study hydrodynamics via AdS/CFT was [1]. The authors related the shear viscosity 
coefficient rj of strongly coupled N = 4 gauge theory plasma in large N limit with the 
absorption cross-section of low energy gravitons by black D3 brane. Other hydrodynamic 
quantities like speed of sound, diffusion constants, drag force on quarks etc can also be 
computed in the context of AdS/CFT (See [2] for review and the references therein). 

If we consider the 't Hooft coupling to be very large but finite, then we have to include 
the string theory contributions to thermodynamic and hydrodynamic quantities, i.e. we need 
to improve the supergravity results by including the higher derivative terms in the action. 
The higher derivative (string theory) corrections to shear viscosity have been calculated in 
[4], [5], [6] 1 . There exists a "viscosity bound conjecture" [3] which states that the viscosity to 
entropy ratio 2 has a lower bound, 

5 > - (1.2) 
s ~ 4vr y 1 

for all relativistic quantum field theories at finite temperature. In fact in presence of the 
R 4 term which is the first higher derivative correction appears in type II string theory the 
viscosity entropy ratio is greater than j-. But the presence of Gauss-Bonnet term in the 
Lagrangian seems to violate the "viscosity bound conjecture" . Interested readers are referred 
to follow [5], [6], [8] for detailed discussions. In all the cases the transport coefficients have 
been determined either by the Kubo formula (graviton absorption) or by quasi-normal mode 
calculation [9], [10] or using the membrane paradigm approach [11]. 

Recently 2 in [12] the authors have developed an elegant systematic framework to con- 
struct the nonlinear fluid dynamics, order by order in boundary derivative expansion. The 
five dimensional Einstein equations with a negative cosmological constant with appropriate 
boundary conditions can be reduced to nonlinear equations of fluid dynamics. In this small 
note we have generalised the construction of local black brane geometry to higher derivative 
gravity. We started with Gauss-Bonnet term as a toy model. We have found the solution 
of a' corrected Einstein equations only up to first order in boundary derivative expansion. 
Once we obtained the a' corrected geometry we calculate the boundary stress tensor up to 
first order in derivative expansion. From the expression of stress tensor one can read the a' 
correction to shear viscosity coefficient. We have also calculated the ratio of shear viscosity 
to entropy and the result agrees with existing results in literature calculated in other ways 
[5], [6]. Though we have found the correction to the metric up to first order in derivative 



^ee also [7]. 

2 Also in [13], authors have discussed second order hydrodynamics for conformal fluid. See also [14], [15] 
and [16] for related discussion. 
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expansion, but it would be interesting to find the corrections to the metric and stress tensor 
up to second order in derivative expansion [27]. 

We proceed in the same way of [12]. The presence of Gauss-Bonnet term preserves all the 
symmetries of AdS^ spacetime. We have also worked in the Eddington-Finkelstein coordinate. 
Our solutions are also non-singular away from r = 0, and specially at the location of the 
horizon. 

The plan of the paper is following. In section 2 we will very briefly sketch the calculation 
frame work of [12] and in section 3 we will display our results. We finish our paper with some 
concluding remarks (section 4). 



2. Fluid Dynamics From Gravity: The Computational Framework 

In this section we will briefly sketch the working procedure of [12]. For detailed discussion 
readers are referred to the original paper. We will also skip the technical details in this 
section. 

• Consider the Einstein-Hilbert action with negative cosmological constant 
where L is the radius of AdS space. 



• The equation of motions are given by 3 



1 6 

Emn = Rmn — -Rqmn — j^Omn = 0- (2-2) 

There exists a class of solutions to these equations of motion given by the "boosted 
black branes" 4 , 

ds 2 = —lu^dx^dr — —^f{br)u ll u v dx^dx v + — ^ {u^u v + rj^) dx^dx u (2.3) 

Li Li 

with, 

u v = -7 

and Ui = 7/3j (2-4) 



where, 7 = 1/y 1 — (5 2 



3 x M = {v,r,x}. 
V 1 = {v,x}. 
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• Putting the values of u^s the metric can also be written as, 

ds 2 = 2jdvdr - —rj 2 j '(br)dv 2 + —rdx l dx l 
L z L z 

+ ^( 7 2 - l)^ 2 - 2 1 f3 l dx l dr - 2^7 2 (1 - f{br))f3 l dx i dv 
+ ^(l-ZW^jW. (2-5) 
The solution is parametrised by four constant parameters b and ft's. 

• The black brane horizon is located at th = 1/6 and the temperature of this black brane 
is given by, 

r =^- < 2 - 6 > 

• Consider the metric 2.5 and replace the constant parameters b and /3j's by slowly varying 
functions 6(x M ) and / Sj(x'*)'s of boundary coordinates 2;^ 

2 2 
(is 2 = 2 7 cfodr - -^7 2 /(&(:r a )r)cfc 2 + ^ H dx i dx i 

2 2 

+ ^ (7 2 - l)dv 2 - 2 1 f3 i {x a )dx i dr - 2-^7 2 (l - f{b{x a )r))f3 l {x a )dx l dv 



r 



2 



+ j^rix - m^ymix^Pjix^dxw. (2.7) 

We will call this metric g i0 \b{x a ) , f3i{x a )) . 

In general the metric 2.7 is not a solution to Einstein equations unless one adds some 
corrections to the metric and also the parameters b(x a ), (3i(x a ) satisfy some set of equa- 
tions, which turn out to be the equations of boundary fluid mechanics. 

Write the parameters b(x a ) and Pi(x a ) and the metric as a derivative expansion of the 
parameters. Up to first order in derivative expansion, 

g = gW (b(x a ), A(z Q )) + (b(x a ),0i(x a )), (2.8) 

b(x a ) = b (0) (x a ) (2.9) 

and 

&(x a )=$\x a ) (2.10) 

where e is a dimension less parameter whose power counts the number of (bound- 
ary)spacetime derivatives acting on the parameters. Since b^\x a ) and [3^' \x a ) do not 
enter in to the first order equation of motions, we have kept the expansion for b and 
fa's up to leading order. 
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• In general one can write the metric and parameters as power series of e. Then plug the 
metric in Einstein equations and solve the metric and the parameters order by order (in 
e). For example in our case since we are interested up to first order, we will plug the 
metric in Einstein equations and solve for and the constraint equations imply some 
relations between the zero th order parameters. We will work in a particular gauge, 

Tr(( 5 (°))-V 1} ) = 0. (2.11) 

• After finding the metric with first order fluctuations one can find the boundary stress 
tensor (using the definition given in [17], [18]). The form of the boundary stress up to 
first order in derivative expansion is given by, 

16irG 5 T^ = T 4 7T 4 L 3 {4u^u u + rj^) - 2T 3 vr 3 L 3 ^, (2.12) 

where a^ u is given by, 

<V = PZPfyw) - \p^d a u a (2.13) 

and P^ v = u^Uy + rj^. 

3. Higher Derivative Correction to First Order Hydrodynamics 

In this section we will explicitly show how one can generalise this procedure to higher deriva- 
tive action. We will consider the Gauss-Bonnet action as a toy example. We will show how 
the first order metric and constraint relations receive a' corrections. As a result the boundary 
stress tensor also receives an a' correction. We have considered the Gauss-Bonnet term as a 
perturbation and hence our metric and stress tensor is correct up to first order in a'. 

3.1 The Action 

We will start with following action, 
where, 

L G b = RmnpqR mnpq - ±RmnR mn + R 2 . (3.2) 
The equation of motion is 

1 6 a' 

EmN = RMN — ^RgMN — -Jj2.9MN — -^9mnLg.b. 

+ 2a' (RmpqlRn QL - 2R pq R M pnq - 2R^R NQ + RRmn) = 0. (3.3) 
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3.2 The Counterterm and Boundary CFT Stress Tensor 

As usual in gravity theories, the action (3.1) should be supplemented with suitable boundary 
terms, for a well-defined variational principle. For Einstein gravity, one consider the Gibbons- 
Hawking surface term [19] 



T E 

ovrG j dM 



l — [ d A x^K , (3.4) 
tG JdM 



where 7^ and K are the induced metric and the trace of the extrinsic curvature of the 
boundary, respectively. 

7n» = 9iiv ~ n^n v , (3.5) 

K^v = -- (V M ra„ + V^) (3.6) 

and is unit outward normal vector to the asymptotic boundary hypersurface. 
A similar term occurs for Gauss-Bonnet gravity and reads [20, 21] 

with E$ is the four-dimensional Einstein tensor of the metric 7^ and J is the trace of 

«V = \ {2KK W KP + K p(T KP°K„ u - 2K^K^K au - K 2 K„ V ) . (3.8) 

Variation of the action I+I b +1^ now gives an expression which does not contain normal 
derivatives of 5g^ u . 

It is well known that the total action has diverges even at tree level. The divergence arise 
from integrating over the infinite volume of spacetime. We regularise this divergence by using 
the procedure proposed by Balasubramanian and Kraus [17], which furnishes a method for 
calculating gravitational action and conserved quantities without reliance on any reference 
spacetime. This technique was inspired by AdS/CFT correspondence and consists of adding 
suitable counter terms I c t to the action of the theory in order to ensure the finiteness of the 
boundary stress tensor [22]. 

We have found that the on-shell action can be regularised by the following counter term 
(see also [23]). 

r 1 /"a,, ,4 ^— ( . , c 2 



8vrG 5 



J dMd 4 x^ (ci + ^ft) (3.9) 



where c\ and C2 are functions of a'. 1Z is the Ricci scalar made out if boundary metric 7. For 
flat boundary geometry (R l x i? 3 ), 1Z = and the counterterm terns out to be (up to first 
order in a'), 



Varying the total action (which contains the boundary terms (3. 4), (3. 7) and (3.10)) with 
respect to the boundary metric 7^, we compute the divergence- free boundary stress-tensor 

<V = - K lw + 2a '( 3 V - J 7^) - J2I1XV + jplit^j ■ ( 3 - n ) 
Since the CFT metric is given by, 

L 2 

= Jim -~r7^> (3-12) 
the boundary CFT stress tensor is given by [18], 



R 2 

T^ u = lim —jS^ 

R^oo *-> 

?2 1 / Q „,/ 



R 1 / 3 a' \ 

= lim JI"^£~ ( K vw ~ K ^ + 2a '( 3 V - J 7^) - + ) ( 3 - 13 ) 

where i? is the cutoff in radial direction. The gauge theory lives on the boundary of the AdS 
space which is at r = R. 

3.3 Solution to the Equation of Motions 

The solution to the equation of motion is given by the metric [24], 

ds 2 = 2dvdr - ^f(br)dv 2 + ^dxW, (3.14) 
where the function f(br) is given by, 

/w=i-(^+^(i+(^). < 3 - i5 > 

The horizon radius th is given by, 

r H = (3.16) 
and the temperature of this black brane is, 

« r »>- 1 'l-'Y (3.17) 



4vr irbL 2 \ L 2 J ' 

By giving a coordinate transformation one can obtain the boosted black brane metric 
which is of the following form {u^ = (—7,7/?)), 

2 9 

y 1 

ds 2 = 2-fdvdr - —r^ 2 f '(br)dv 2 + —,dx % dx % 

2 2 
+ ^( 7 2 - l)dv 2 - 2 1 [3 i dx i dr - 2^ 7 2 (1 - f (br)) fadx* dv 



r 



2 



+ - /(6r))ft/^xW. (3.18) 
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The CFT metric is given by 3.12. In a' — > limit the CFT metric h^ v is simply the 
Minkowski metric r)^ v . But when we include the a' correction then the CFT metric is no 
longer 77^. It has the following form (for the metric 3.14), 

o\ 
10 

10 

1/ 

For metric 3.18, h^ v also has off-diagonal terms proportional to a'. If we want to keep the 
CFT metric to be rjn U then we have to rescale either the time coordinate v — > (1 — a'/L 2 ) v 
or the space coordinates x l — > (1 + a' /L 2 ) x l . Here we will rescale the time coordinate 



\ 



(3.19) 



In this rescaled coordinate (V,x) the metric 3.14 takes the following form, 



ds 2 



2 ( 1 ~ Jl ) dVdr ~ j2^ hr ) dV + j^dx*dx % . 



And the temperature has also been rescaled to, 

1 



irbL 2 



1 



2a' 

T 2 



(3.20) 



(3.21) 



(3.22) 



Now we can give a coordinate transformation (boost) and the boosted black brane metric in 
(V,x) coordinate becomes, 

ds 2 



2- ( 1 - ^) dVdr - ^ 2 f(br)dV 2 + ^dx l dx l 



+ — h 2 



l)dV 2 - 2 7 1 



J Pidx'dr - 2-^7 (1 " f(br))(3 t dx l dV 



where 



+ 1 2 j^(l-f(br))p i p j dx i dx j 
f(br) : 



(3.23) 
(3.24) 



3.4 First Order Metric 



In this section we will determine the metric to the first order in the derivative expansion. The 
metric and the parameters are given by Eq. 2.8, 2.9 and 2.10. 

We will choose a coordinate to set = { — 1,0,0,0} (i.e. uy = —1, (3^ = 0,7 = 1) and 
6(°) = 1 at a given point x^ = . Around this point the velocities and the temperature fields 
can be expanded up to first order in derivative, 



3 From now, our = {V, x). 
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6 = 1 + x^6 (0) (3.25) 

and 

$ 0) = x»d^f ) . (3.26) 
So the metric up to first order in derivatives is given by, 

/ / \ 2 2 

ds 2 = 2 ( 1 - ^ ) dVdr - ^f{r)dV 2 + ^dx i dx i 



2 { l ~ i) x ^ d ^t )dxldr 



2 3 

- 2-^(1 - f(r))x»d^f ] dx l dV - ^d^f{r)dV 2 . (3.27) 

This is the metric g^ expanded up to first order in derivative. 

As we explained earlier is not a solution to equation of motion. We have to find the 
fluctuation metric g^> , which added to </°) solves the equation of motion up to first order in 
derivative expansion. 

Because of spatial SO (3) symmetry of background black brane metric, we can separately 
solve for the SO(3) scalars, SO(3) vectors and SO(3) symmetric traceless components of 

3.4.1 The Scalar Sector 

The scalar components of gW are parametrised in the following way, 

9u ( r ) = ^T2^i( r ) {sum over i) 
I? 

9vU r ) = 72- fc i( r ) 

9 { vliT) = -l{l-^)h l {r). (3.28) 

Here we are working in the following gauge, 

Tr(( 5 (°))-V 1} ) = 0. (3.29) 

The scalar Einstein equations (equations invariant under SO(3) rotations) are divided up 
into constraint and dynamical equations. 

Constraint Equation 1 

The first scalar constraint is 

2 

^f(r)E Vr + E vv = 0, (3.30) 
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which evaluates to 



(l-|Q(^<»>-i%f>)=0 



i.e. d v b< a > -jd,^ =0. (3.31) 

The constraint relation remains unchanged in presence of Gauss-Bonnet correction. This 
relation is a consequence of the conservation of boundary energy momentum tensor. 

Constraint Equation 2 

The second constraint equation is 

j^f{r)E rr + E Vr = 0, (3.32) 



leads to 



127-% (r) + (3r 4 - l)h[ (r) - L% (r) = -2L V ( 1 - ) $/?< 0) . (3.33) 



Dynamical Scalar Equation 

In addition to these constraint equations we have to add one dynamical equation. Like [12] 
we will add the following simplest equation, 

E„ = 5h[ (r) + rh'[(r) = 0. (3.34) 

We will solve the equations 3.33 and 3.34 to find out the function h\(r) and k\(r). These 
two equations satisfy all the regularity and normalisation conditions explained in [12]. Hence 
the solutions are given by, 



/ii(r) = 
2 

3L 2 L\ 

The constants appear in h±(r) and k±(r) can be set to zero following the same argument in 
[12]. So the scalar part of the fluctuation metric is given by 6 , 

(gW)^ dx«dx? = \r(±- fi) Wf* dV\ (3.36) 



e '9s>9vi9T are the scalar, vector and tensor part of the fluctuation metric g 1 respectively. 
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3.4.2 The Vector Sector 

We will parametrise the vector part of the fluctuation metric as, 

Constraint Equation 3 



($) = 2^(1 - f~(r))j^dVdx\ (3.37) 



In this sector the constraint equation is given by, 

„2 



^f(r)E ri + E Vi = (3.38) 
which gives, 

(*6»>-M») = 
i.e. dM 0) - c\f3f ] = 0. (3.39) 

Again the presence of Gauss-Bonnet term does not have any effect on this constraint equation 
which can be interpreted as a consequence of conservation of boundary stress tensor. 

Dynamical Equation for l-^(r) 

The dynamical equation in vector sector is, E„ = 0. The equation for ][\r) turns out 
to be 7 , 

dr ^r 3 dr J% y 1 ) r 2 v M i jji \^ r 6 r 2 

The solution is given by, 



i ( -Aj"(r)) = -3^M 0) + £ ( ^ " "4 ) W?\ (3.40) 



j| 1) (r) = LV 



L 2 \ r 4 



^/3f 0) . (3.41) 



The integration constants are set to zero since the solution should be normalisable at the 
boundary and the stress tensor must be renormalisable [12]. So the vector part of the fluctu- 
ation metric is given by, 

(s^) ^ dx a dx p = 2rd v l3f ] (l - |Q dVdx\ (3.42) 



7 when j^'s or their derivatives appear with a' then we have set the leading order (in a') values of j^'s 
or their derivatives. 
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3.4.3 The Tensor Sector 

The SO(3) tensor part of can be parametrise in the following way 

' a/3 



(#) = pogWcfc* (3.43) 



where af^ is symmetric traceless 3x3 matrix. 



Dynamical Equation for a[j\r) 



Equation for (r) follows from the Eij = 0. 

d ( ( 5 \ d \\ rt2 2 (0) 

fo[(r -r)—a\ j '(r)j=-6Lr < 7^ 

(2(9 + 5r 8 )«;(. 1) (r)) 

L 2 r 4 



+ ^4 (2L 2 r 2 (-4 + 9r 4 )4 0) + 2r(-3 + r 8 )a"/ 1} (r))(3.44) 



where, 

^=¥?lM^f (3-45) 
The solution for off (r) is given by 8 , 

- 2 (r (l - - ^ (l - ff)) (346) 

Summary of a' corrected first order calculation 

The a' corrected metric </°) + expanded up to first order in boundary derivatives 
around some point where = 1 and = { — 1,0,0,0} is given by, 

/ / \ 2 2 

ds 2 = 2 ( 1 - %r ] dVdr - ^-rf(r)dV 2 + ^dx i dx i 



2 ( 1 ~T 2 ) x ^ d ^f )dxldr 



2 3 

2-^(1 - f(r))x»d^f ] dx l dV - y x ^b^~f{r)dV 2 



+ \r(l- ^) d4 0) dV 2 + 2rdy^ (l - ^) dVdx* 

^M'-vYM 1 -^)**^- (3 - 47) 



3 In the right hand side of Eq. 3.44 we set the leading order (in a') values of a'j and their derivatives. 
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Global solution to first order in derivative 



The metric 3.47 has been calculated about x M = assuming &(°) = 1 and [if ] = 0. But 
one can also write the metric about any point. The a! corrected global metric is given by, 

(Cx! \ T 2 T 2 

1 — -j-p, J u^dx^dr — -j-pj '(br^^Uvdx^ l dx u + j-^P^dx^dx" 

+ 2r 2 bF(br)cr^dx^dx v + - 1 - d x u x dx^dx v 

~ r i 1 ~ I 1 ) ^M^u^dx^dx" (3.48) 

where F{r) is given by, 

*M=(i(i-£)-^(i-£)). 

3.5 The Energy Momentum Tensor 

Once we obtained the 0/ corrected metric up to first order fluctuations, we can find the 
boundary stress tensor using the definition 3.13. Different components are given by, 

16nG 5 T Vi = -± f 1 _ |Q x ^ 0) 

-4TWfl + ^x^/3f. (3.50) 



L 2 , 

Here we have written the energy momentum tensor about for the metric 3.47 i.e. about 
x^ = assuming that = 1 and (3^ = at the origin. But for the global metric 3.48 the 
energy momentum tensor can be written in a covariant way. Up to first order in derivative 
expansion the energy momentum tensor has the form given by Eq. 1.1 (see [26]). Therefore in 
presence of the Gauss-Bonnet term one can write the energy momentum tensor in a covariant 
form in the following way, 



11a' 

IT 1 



= T\ 4 L 3 (l + |Q (4u,u u + 7^) - 2TVL 3 (l - |Q a^. (3.51) 
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The shear viscosity coefficient is given by, 



' 16ttG 5 L 3 6 3 V 
Entropy density is given by, 

S 



where V3 = f d 3 x. Hence, 

77 

s 

which is in agreement with [5], [6]. 
4. Discussion 

In this paper we have constructed the local black D3 brane geometry in presence of Gauss- 
Bonnet term in the bulk action up to first order in derivative expansion. The local solution we 
found is non-singular except at r = 0. We used the counter term method of Balasubramanian 
and Kraus to find the boundary stress tensor , which we have expanded up to first order in 
derivatives, for the a' corrected metric. From the expression of stress tensor one can read the 
shear viscosity coefficient. We have found the a' correction to the shear viscosity coefficient 
and also to the viscosity entropy ratio. These results are in agreement with the existing 
results in the literature. 

Recently in [25], authors have demonstrated that given a black brane geometry with 
regular event horizon the location of the horizon in radial direction turns out to be a local 
function of fluid dynamical variables evaluated at the corresponding points on the boundary. 
In presence of this regular local event horizon they constructed an appropriate area form on 
spatial section and then taking the pull-back of this area form to the boundary they defined 
a local entropy current for the dual field theory. The entropy current has the following form 

4G 5 L 3 & 3 j£ = u» + 0{e 2 ). (4.1) 

The first order (0(e)) correction to entropy current is zero. As mentioned in [25] finding 
out the a' correction to entropy current would be an interesting problem to solve. Since the 
relations Eq. 3.31 and Eq. 3.39, which follow from the conservation of energy momentum 



L 2 J 16ttG 5 V L2 



Area 



4G 5 V 3 



4L 3 G, 



5 

I ( x _ 3a' 



4vr V L 2 



4L 3 6 3 G 5 V L 2 

TW 'l + £ I • (3,3) 



4G 5 V L 2 
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tensor, remain unchanged (up to an overall factor), therefore, given the a' corrected bulk 
metric 3.48 it is easy to check that the entropy current, up to 0(e), has the following form , 



Again the entropy current does not receive any 0(e) correction. So to find out the a' correction 
up to order 0(e 2 ) to the entropy current, one has to first find the corrected metric up to 0(e 2 ) 



Although we have considered only the four derivative terms in the Lagrangian but it 
would be very interesting to generalise this idea for any higher derivative gravity. From the 
string theory point of view it would be nice to construct the corrected local black D3 brane 
geometry for R 4 term [28]. 
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